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^ , Abstract 

^ . It is known that random convex polygonal lines on (with the endpoints fixed at = 

I (0, 0) and n = (ni, 712) — )■ 00) have a limit shape with respect to the uniform probability 

Q ' measure, identified as the pai^abola arc ^ c(l — 2:1) + ^Jx2 = \fc, where n2/ni — )■ c. The 

present paper is concerned with the inverse problem of the limit shape. We show that for 
^ \ any strictly convex, C^-smooth arc 7 C starting at the origin, there is a probability 

measure Pji on convex polygonal lines, under which the curve 7 is their limit shape. 
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1. Introduction 



Consider a convex lattice polygonal line F with vertices on := & 1^ ■ i, j > 0}, 

• starting at the origin and such that the slope of each of its edges is nonnegative and does 

not exceed the angle of 90°. Convexity means that the slope of consecutive edges is strictly 
increasing. Let 77 be the set of all such polygonal lines with finitely many edges, and by i7„ 
the subset of polygonal lines F E n with the right endpoint fixed at n = (ni, ^2) G Z^. 

The limit shape, with respect to a probability measure P„ on i7„ as n -> 00, is understood 
as a planar curve 7* such that, for any e > 0, 

b : 1™ Pn{r e i7„ : c?(r„,7*) < 4 = i, (i.i) 



where Fn = Sn{F), subject to a suitable scaling S'„ : — )■ M^, and d(-, ■) is some metric on 
the path space, e.g., induced by the Hausdorff distance between compact sets, 



d-^iA, B) := max < maxmin |x — yl, maxmin |x — y| > . (1-2) 

xeA yeB y£B xgA J 

Of course, the limit shape and its very existence may depend on the probability law P„. 
With respect to the uniform distribution on i7„, the problem was solved independently by 
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Vershik [[T4l . Barany |l2l and Sinai [|T2ll . who showed that if rii, n2 — c>o so that n2/ni c G 
(0, oo), then under the scaling S'„ : (xi, X2) ^ {xi/ni,X2/n2) limit (11.11) holds with respect 
to the Hausdorff metric dn and with the limit shape 7* identified as a parabola arc 

c(l - Xi) + ^/X^ = y/c, 0<Xi,X2<l. (1.3) 

Recently, Bogachev and Zarbaliev |I71[8]| proved that the same limit shape (11.31) appears for 
a large class of measures _P„ of the form 

P^in := H h^, B^:= U K ^ il.), (1.4) 

where the product is taken over all edges Cj of -T G iln, /ci is the number of lattice points on 
the edge Cj except its left endpoint, and 

/r + k-l\ r(r + 1) ■ ■ ■ (r + k - 1) 

, )=J 4^ K k = 0,l,2,.... (1.5) 

This result has provided first evidence in support of a conjecture on the limit shape universal- 
ity, put forward independently by Vershik lfT4l p. 20] and Prokhorov [[TTI . The class of prob- 
ability measures (11.41) with coefficients (11.51) belongs to a general meta-type of decomposable 
combinatorial structures known as multisets (see [[B §2.2]). Bogachev H has extended the 
universality result to a much wider class of multiplicative probability measures (11.41) including 
the analogues of two other well-known meta-types of decomposable structures — selections 
and assemblies (cf. [1, §2.2]); for example, this class includes the uniform distribution on the 
subset of "simple" polygonal lines (i.e., those with no lattice points apart from vertices). 

However, universality of the limit shape 7* given by (11.31) has its boundaries; indeed, in the 
present paper we consider the inverse problem and show that any -smooth, strictly convex 
arc 7 G (started at the origin) may appear as the limit shape with respect to a suitable 
probability measure on i7„, as n — )■ 00. For early drafts of this result (treated in terms of 
approximation of convex curves by random polygonal lines), see 

Like in dUQ, our construction employs an elegant probabilistic approach based on ran- 
domization and conditioning (see Ol), first used in the polygonal context by Sinai p2l. The 
idea is to represent a given measure P„ on 7J„ as the conditional distribution, Pn(P) = 
Q{r I Un), induced by a suitable "global" probability measure Q defined on the space U = 
U„i7„ of convex lattice polygonal lines with a free right endpoint. In turn, the measure Q = Qz 
depending on a two-dimensional parameter z = {zi, Z2) is constructed as the distribution of a 
suitable integer- valued random field v = z/(-) with mutually independent components, defined 
on the subset A" C consisting of points X — (3^!) X2 ) G with co-prime coordinates. 
Note that a polygonal line P G i7 is easily retrieved from a configuration {i^{x)}xex using the 
collection of the corresponding edges {xu{x)}xex and the convexity property. 

It turns out, however, that in order to fit a given curve 7 the parameter z = (zi, Z2) needs 
to allow for a dependence on a; G X. We derive suitable parameter functions zi{x) and Z2{x), 
assuming that they depend on x through the ratio X2/X1 only, which is particularly convenient 
in conjunction with the parameterization of the curve 7 using its tangent slope. As one would 
anticipate, if 7 = 7* then the functions zi{x), Z2{x) are reduced to constants and our method 
recovers the uniform distribution on 

To summarize, our main result is as follows. 
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Theorem 1.1. Let 7 C be a strictly convex, -smooth arc, with the endpoints and (1, c-y) 
and with the curvature bounded from below by a positive constant. Suppose that n^jrix — )■ c^, 
and set Fn := n'[^r. Then there is a probability measure on i7„ such that, for any e > 0, 

lim P2{r e Un : dniPn, l)<e] = 1. (1.6) 

n^oo 

Remark 1.1. It will be more convenient to use another metric on the space of convex paths 
(denoted by dc), based on the tangential parameterization of paths and a sup-distance between 
the corresponding arc lengths. However, the metrics dc and d-^ are equivalent. 

Remark 1.2. It is interesting to try and relax the C'^-smoothness condition on 7 (e.g., by per- 
mitting "change-points" or corners), as well as to allow for degeneration of the curvature (e.g., 
through possible flat segments). We will address these issues elsewhere. 

The rest of the paper is organized as follows. In Section [2l we introduce the space of 
convex paths on the plane and endow it with a suitable metric. In Section [3l the measures Q1 
and P2 are constructed for a given convex curve 7. In Section HI the parameter vector-function 
z{x) is chosen to guarantee the convergence of "expected" scaled polygonal lines Pn = n^^P 
to the target curve 7 (Theorem 14.21) . Refined first-order moment asymptotics are obtained in 
Section [51 while higher-order moment sums are analyzed in Section [6l Section |7] is devoted to 
the proof of a local central limit theorem (Theorem 17.11) . Finally, the limit shape result, with 
respect to both Q1 and P^, is proved in Section [8] (Theorems 18.11 and 18 . 21 respectively). 

Some notations. For a row- vector a; = (xi, X2) G M^, its Euclidean norm (length) is denoted 

by |x| := {x\ + xlY/"^, and (x, y) := xy^ = xiyi + X2?/2 is the corresponding inner product 
of vectors x,y e M^. We denote := {k e Z : k > 0}, := Z+ x and similarly 
M+ := {x e M : a: > 0}, := M+ x M+. For 2 = (21,22) e andx = (xi,X2) G Z^, we 
use the multi-index notation := Zi^Z2'^. The notation a„ x 6„ as n — )■ 00 means that 

< liminf — < lim sup — < 00. 

yn n^oo y-n 

We also use the standard notation x„ ~ y„ for Xn/Vn — > 1 as rii, 77-2 — 00. 



2. Preliminaries: convex planar paths 

Let (yf be a bounded function defined on some interval [0, a], such that (7(0) = 0, and suppose 
that g is non-decreasing and convex on [0, a]. Convexity means that the function's epigraph 
{(u, w) : < n < a, g{u) < t;} is a convex set on the plane. Furthermore, assume that g is 
continuous on [0, a] and piecewise differentiable, with the derivative g' continuous everywhere 
except at a finite set of points (we allow g\a) to be infinite, g'{a) < +00). It follows that 
the function t = g'{u) is nonnegative and non-decreasing in its domain, and in particular 
< to < 9'{^) ^ ^1 ^ 00 (0 < n < a), where 

to := inf g'{u), h := sup g'{u). (2.1) 

0<u<a 0<u<a 

Denote by 7^ = 7 the graph of a function g with the above properties, and let (5 be the set 
of all such curves. For the spaces i7„, 77 of convex polygonal lines introduced ablve, we have 
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Un C n C 6. If a polygonal line F E Iln is taken as a "curve" 7, then the corresponding 
function g = gr is a piecewise linear function defined on [0, ni], such that gr{ni) = n2. 

Let us now equip the space 25 with a suitable metric. If a function g = g-y determines a 
convex curve 7 G 35, we set 

u^{t) := snp{u : g'^{u) <t}, < t < 00, (2.2) 

with the convention that sup = 0. That is, u^(t) is a generalized inverse of the derivative t = 
g'^{u) (cf. llll §1.5]). It follows that the function u^{-) is non-decreasing and right-continuous 
on [0, 00], with values in [0, a]. Moreover, if to, ^1 are the extreme values of the derivative g'^ 
(see (12.11) ) then u^{t) = for all t < and Uj{t) = a for all t >ti. 

Denote by f'y(t) the length of the part of 7 where the tangent slope does not exceed t, 

u~, it) 



{t)= Jl + g'Juydu, 0<t<oo. (2.3) 



Our assumptions imply that every curve 7 G is rectifiable, since 



^7V 



00) = y 1 + g'^{uy du < J {1 + g'^{u))du = a + g^{a) < 00. 

Finally, we define the function dc ■ & x & ^ IR+ U {+00} by setting 

c?£(7i,72):= sup Kit) - i^,{t)\, 71,72 G 25. (2.4) 

0<t<oo 

Proposition 2.1. The function dc{-, ■) satisfies all properties of a distance. 

Proof Clearly, ^£(71, 72) = ^£(72, 7i) and ^£(7, 7) = 0. The triangle axiom is also obvious. 
So it remains to verify that if ^£(71,72) = then 71 = 72. To this end, approximating 
7i, 72 G 25 by C^- smooth strictly convex curves 7^ , 7I, respectively, we reduce the problem 
to checking that if 7^ , 7I are close to each other in the sense of dc, then they are also close in 
the Hausdorff metric rf^; that is, if dd-^i, 72) — ^ then rf^(7i , 73) — (A; — )■ 00). 

Next, for a strictly convex, increasing function g^ G C^[0, a], the function u^{t) defined in 
(12. 2|) is given explicitly by 

r 0, < t < to, 

«7W = \ to<t<h, (2.5) 

y a, ti < t < 00, 

where {g'^)~^{t) is the (ordinary) inverse of the derivative g'^{u). In particular, the equations 
u = Uj{t), V = g^{u^{t)) determine a parameterization of the curve 7 via the derivative 
t = g!y{u). Differentiating formula (12.31) with respect to t, we find 



du.y 1 d£j 

'df ~ yrTt2 It" 



u^(0) = 0, (2.6) 



and hence 



^ = ^.^ = ^^^, .,(0) = 0. (2.7) 
dt du dt y/TT^ dt ' ^ 
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Integrating equations (12.61) . (12.71) by parts, we obtain 



s 



Note that these equations are linear in Hence, setting for 7^^, 7I 

AMfc(t) := u^k{t) - u^k{t), Avk(t) := v^k{t) - v^k{t), 
A4(t) :=i,k{t)-i^k{t), 

from (|2.8|) we get 

This implies that if Aikit) — > 0, uniformly in t G [0, 00], then Auk{t) — )■ 0, Ai;fe(t) — )■ 0, also 
uniformly on [0, 00] (k — )■ 00). This completes the proof. □ 
From the proof of Proposition [27T1 one can see that the following result holds. 

Corollary 2.2. The metrics dc and d-^ are equivalent; in particular, dci'jk, 7) — ^ if and only 
if dnilk,l) 0. 

Consider a fixed convex curve 7 G ©, represented as the graph of an increasing convex 
function g^, which for definiteness is assumed to be defined on the interval [0, 1]. We will be 
working under the following 

Assumption 2.1. The function is strictly increasing and strictly convex on [0, 1], and g^^ E 
C^[0, 1]. In particular, g'yiu) > 0, g"{u) > for all u e [0, 1]. Moreover, the curvature of 
the curve 7, given by the formula 



9» 

(1 + 

is uniformly bounded from below. 



X7H = 7T— TT^wI' 0<n<l, (2.9) 
(1 + ^;(n)2)3/2 



inf kJu) > Ko > 0. (2.10) 



As was mentioned in the proof of Proposition 12. 1[ the graph 7 of the function g^ can be 
parameterized by the derivative t = g!y{u) via the equations u = u^{t), v = g^{u^{t)), where 
u^{t) is given by (12.51) . Expression (12.91) for the curvature is then reduced to 

g''{uJt)) 

><7(^) = (7^;,y;, , ^o<t<ti, (2.11) 

whereto = inf^gf' (m), U = sup^ g'Ju) (see (O))- 



5 



3. Construction of the measures Q], and 



Consider the set A" C of all pairs of co-prime nonnegative integers, 

X := {x= (xi,X2) e : gcd(xi,X2) = 1}, (3.1) 

where "gcd" stands for "greatest common divisor". Denote by r(x) := X2/X1 E [0, +00] 
the slope of the vector x = {xi,X2) E X. Let <P := (Z^)"^ be the space of functions on X 
with nonnegative integer values, and consider the subspace of functions with finite support, 
'pQ := {u E ^ : #(supp i^) < 00}, where supp i' := {x E X : i^(x) > 0}. It is easy 
to see that the space <Po is in one-to-one correspondence with the space IJ = Unez^ -^n of 
all (finite) convex lattice polygonal lines, whereby each x E X determines the direction of a 
potential edge, only utilized if x G supp u, in which case the value z/(x) > specifies the 
scaling factor, altogether yielding a vector edge xu{x); finally, assembling all such edges into 
a polygonal line is uniquely determined by the fixation of the starting point (at the origin) and 
the convexity property. Note that ^{x) = formally corresponds to the "trivial" polygonal 
line with coinciding endpoints. In what follows, we identify the spaces 77 and ^q{X). 
Let us now introduce on '^o{X) a probability measure by setting 

QUn := n ^"'^"Hl - ^1, UbT^uE M^), (3.2) 

where z = z{x) = {zi(x), Z2{x)) is a parameter (vector- )f unction such that < Zj{x) < 1 
(x E X). Its explicit form, determined by a given curve 7 G (S, will be specified later on. So 
far, we only assume that 

Z:= l[{l-z^)>0, (3.3) 

which guarantees that (13.21) is well defined. Definition (13.21) implies that the random variables 
{i'{x)}xex are mutually independent and have geometric distribution with parameter z^, 

Q2{u{x) = k} = z'^^'il - ^"), k E Z+; (3.4) 

in particular, the corresponding expected value and variance are given by BH §XI.2, p. 269] 

ElHx)] = = E Var[.(x)] = = kz'^. (3.5) 

k=l ^ ' k=l 

Note that Q1 can be extended in a standard way to a measure on the space '^{X) of all 
nonnegative integer- valued functions on X. However, Q1 is in fact concentrated on the subset 
C ^(A") consisting of all finite configurations z/(-). 

Lemma 3.1. Condition (13.31) is necessary and sufficient in order that Q1{i' E 'Pq{X)} = 1. 

Proof. According to (13.41) . 

J2 q:m^) > 0} = < 00 

x&x xex 

whenever the infinite product in (13.31) is convergent. By Borel-Cantelli's lemma, this implies 
that only finitely many events {i^(x) > 0} may occur (QJ-a.s), and the lemma is proved. □ 
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As a result, with -probability 1 a realization of the random field determines a (ran- 
dom) convex polygonal line F E 11. Denote by = {^i, C2) the right endpoint of F, so that 
Un = {F E n : = n}. Accordingly, QJ induces the conditional distribution on i7„, 

O'^(F) 

P'nin := Ql{F I ir = n} = ^.Tf^J . , T G i7„. (3.6) 



4. The choice of the parameter function z{x) 

In the above construction, the measure P] depends on the vector parameter {z{x)}xex- So far, 
this function was only assumed to guarantee the convergence of the infinite product in (13.31) . 
Let us now adjust it to a given curve 76©. 

Let F{t) denote the part of the polygonal line F E U where the slope of edges does not 
exceed t E [0,oo]. Set X{t) := {x E X : t{x) < t}. Recalling the association F -h^ u 
described in Section [3l the polygonal line F{t) is determined by the truncated configuration 
v{x)lx{t){x). Denote by ^(t) = (^i(t), ^2(i)) the right endpoint of r(t), 

e,(t)= ^:'-^(^) 0' = 1'2), (4.1) 

x&X{t) 

and by ^r{t) its length, 

ij.{t) = J2 (4.2) 

Let us impose the following calibration condition, 

lim n:[^E][£r{t)] = iM), 0<t<oo, (4.3) 

n— >oo 

where E] stands for the expectation with respect to the measure QJ and iy{t) is the corre- 
sponding length function associated with a given curve 7. More specifically, denote p„ := 

Cy/cn 1 and set x := (xi, p„X2) for x = (xi, X2). We will seek zi{x), Z2{x) in the form 

Zj{x) = exp {-anSjirii))} (j = 1,2), (4.4) 

where 

an := {Puni)'^^'^ 0, r(x) = X2/X1 = p„X2/xi, (4.5) 
and 6{t) = (5i(t), 52{t)) is a function on [0, 00] such that 

inf 6 At) >6,>0 (j = l,2). (4.6) 

0<t<oo 

Remark 4.1. Note that the right endpoint of the scaled polygonal line Fn := n^^F {F E IT) 
has the coordinates (1, c„), where c„ := n2/ni, whereas the right endpoint of the arc 7 lies at 
the point (1, c^), where := (?-y(l) (0 < < 00). Hence, in order for relation (11.61) to be 
true, it is natural to pass to the limit n — )■ 00 in such a way that c„ — )■ c^. In what follows, we 
will always be assuming that this condition is fulfilled. 
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According to (13.51) and (14.21) . and using notation (14.41) . we have 

oo oo 

El[^r{t)]= I^E^'" = 5Z Yl Ixle-""'^^^'^^"^*))^ (4.7) 

xGX{t) k=l k=l x£X(t) 

To deal with sums over the sets X(t) C X, the following lemma will be instrumental. 
Recall that the Mobius function /i(m) (m G N) is defined as follows: /i(l) := 1, /i(m) : = 
(— l)*^ if m is a product of d different prime numbers, and /i(m) := otherwise (see IfTOl 
§16.3, p. 234]); in particular, |/i(m)| < 1 for all m G N. 

Lemma 4.1. Let / : — )■ M a function such that /(O, 0) = and 

oo 

J] J] |/(/iA;x)| <oo, h>0. (4.8) 

k=l x&Zl 

For h > 0, consider the functions 

oo 

m=l x-gA' xGA' 

r/zen the following identities hold for all h > 

oo 

F{h) = J2 /(M, ^H^) = (4.10) 

Proof. Recalling definition (13.11) of the set X, observe that = |Jm=o "^'^^ hence, definition 
(1431) of F(-) is reduced to (14.101) . Representation (14.101) for F'^{-) follows from the Mobius 
inversion formula (see [10, Theorem 270, p. 237]), provided that J2k m \F\hkm)\ < oo. To 
verify the last condition, using (14.91) we obtain 



oo / oo 



\F^{kmh)\ <Y[J2Y1 \f{hkmx)\ = J] E < 

k,m=l k=l \m=lxGA' / k=l xeZ^ 

according to □ 

Theorem 4.2. Suppose that the functions Si{t), 62{t) satisfy condition (14.61) . Then, in order 
that equation (14.31) be fulfilled for all t G [0, oo], it is necessary and sufficient that 

5,(t) = +oo (j = l,2), t<to, t>ti, (4.11) 

5i{t) + t52{t) = Kg'!^{u^{t)f'\ to<t< ti, (4.12) 

where k := (2C(3)/C(2))^/'^, C(s) := Yl'kLi '^/^^ the Riemann zeta function, and the func- 
tion u^(t) is given by (12.51) . 

Proof. Let us set 

fix) := |x| e-""<*'^«*») l;t(i)(a;), x G M^, (4.13) 



8 



suppressing for simplicity the dependence on t and n. Following notations (14.91 ) of Lemma 
14. 1[ equation (|4.7I) is rewritten in the form 



E:[irit)]=J2k~'F\k), (4.14) 



k=l 



whereas from (14.101) we have 



F{h) = J2 Yl hlxle-""^^^'^^^^^^^^ (4.15) 

a;i=0 0<X2<txi 

oo oo 

xi,X2=0 y=0 

^ 0{l)a-^h-^. (4.16) 
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In particular, this gives F{hk) = 0{k"'^), uniformly in A; G N, and it follows that condition 
(14.81) of Lemma l4.1l is satisfied. Hence, using (14.101) and (14.151) and recalling that n^^ = Pnd^, 
from (|4.14|) with h = k we obtain 



Y^E2[£r{t)]= Pnal J] mfi{m)F{km) 

k,m=l 
00 00 

= Pn '^f'^'^^Yl Y an|2;|e-'=™""<"''^("("»^ (4.17) 



k,m=l xi=l 0<X2<txi 



Taking into account estimate (|4.16l) . we see that the general term in the double sum over 
k, m in (14.171) admits a uniform bound of the form 0{1) k'^m'"^, which is a term of a conver- 
gent series. Therefore, we can apply Lebesgue's dominated convergence theorem to pass to 
the limit in (14.171) termwise, as n — )■ 00 (i.e., — )■ 0). In order to find this limit, note that the 
internal double series over xi, X2 in (14.171) is a Riemann sum for the integral 



Q<.X2 <tx-i 



yjx\ + xl e-'=-(-i'5i(-2/-i)+-252(x2Ai)) (418) 



Moreover, this sum does converge to integral (14.181) as a„ — 0, since the integrand function 
in (14.181) is directly Riemann integrable, as follows from an estimation similar to (14.161) . 
By the change of variables y\ = u,y2 = us integral (14.181) is reduced to 



ft / POO 



kmuiS,is)+sS2is)) d« ) = / + ^ ds. 



{kmy 7o (5i(s) + s62{s)y 



Substituting this into (14.171) we get 



lim nT^E]Mt)] = 2 E A E ^ 3 



'^^'^ ^ d., (4.19) 



C(2) Jo {6i{s) + s62{s)y 
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where we used the identity Ylm=i ^l^ifn) = C(2) ^, which readily follows by the Mobius 
inversion formula dHOl) applied to F^{h) = h~^, F{h) = Em=i(^^)"^ = ^^"^((2) (cf. EOl 
§17.5, Theorem 287, p. 250]). Recalling the notation k introduced in Theorem 14.21 and using 
condition (14.31) . from (14.191) we obtain 



/ 3 ds = L,{t) , < t < cx). (4.20) 

Jo {Si{s) + s52{s)) 

According to definitions (12.21) and (12.31) . we have ij{t) = for t G [0, to) and i^{t) = £^(oo) 
for t G (ti,oo], while for t G (to^^i) the derivative i'^{t) is determined by formula (12.61) . 
Hence, differentiating identity (14.201) with respect to t, we obtain (14.1 II) and (14.121) . □ 
Let us now check that equation (14.121) has a suitable solution. 

Proposition 4.3. For t G [to, ^i] l^t "■^ 



5,{t) := «:x,(t)V3 ^1^^^, S,{t) := (4.21) 

where = g^{l) and the curvature x^it) is given by (12.1 II) . Then the functions Si{t), 62(t) 
satisfy assumption (14.61) and equation (14.121) . 

Proof. It is straightforward to verify that equation (14.121) is satisfied. A lower bound of the 
form (14.61) follows from assumption (12.101) . □ 

Remark 4.2. In the "classical" case, where the curve 7 = 7* is determined by equation (11.31) . 
it is easy to check that the corresponding curvature (see (12.91 )) is given by 

cil + t/cf 

Hence, expressions (14.211) are reduced to the constants 5i = k(c/2)^/3, 82 = 5i/c (cf. f/l). 

Assumption 4.1. Throughout the rest of the paper, we assume that the parameters zi{x), Z2{x) 
{x G X) are chosen according to formulas (14.41) with the functions 5i{t), 52{t) given by (14.1 II) . 
(14.211) . In particular, the measure QJ becomes dependent on n = (721,722), as well as the 
-probabilities and the corresponding expected values. 



5. Asymptotics of the expectation 

In this section, we derive a few corollaries from the above choice of zi{x), Z2{x), assuming 
throughout that Assumptions 12. II and 14. II are satisfied. 

Theorem 5.1. The convergence in (14.31) is uniform in t G [0, 00], 

lim sup |n^^E2'[£r(t)] -^7(t)| = 0. (5.1) 

n-)-ooo<t<oo 

We will use the following simple criterion (see |l7l Lemma 4.3]). 

Lemma 5.2. Let {/„(t)} be a sequence of non-decreasing functions on a finite interval [a, h], 
such that, for each t G [a, 6], lim„_j,oo /n(t) = fit), where fit) is a continuous {non- 
decreasing) function on [a, 6]. Then fn{t) — > /(t) uniformly on [a,b]. 
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Proof of Theorem [STTI Note that for each n the function 

Ut) := ni'E][irit)] = - l^l^lH^)] 

rii ^ — ' 

is non-decreasing in t and the limiting function f(t) := i^{t) given by (I2.3h is continuous on 
[0, oo]. Hence, by Lemma [5^ the convergence in (15.11) is uniform in t on every finite interval 
[0, t*]. To complete the proof, it suffices to check that for any e > and for large enough n, 
there exists t* < oo such that for all t >t* 

ni^E][£r{oc) - ir{t)] < e. (5.2) 

Using (14.71) . similarly to (14.161) we can write 

oo 

k=l x&X\X(t) 
oo oo 

^ (xi+X2)e-""*^'^*("^+"^)/^ (5.3) 

fe = l X\=l X2>tXl 

Note that the number of integer pairs {xi,X2) (with xi > 1, X2 > 0) satisfying the conditions 
xi + X2 = y and X2 > txi does not exceed y/{t + 1). Hence, again using estimate (14.161) . we 
see that the right-hand side of (15.31) is bounded from above by 



oo oo 



V V ^ = ^ V ^ = ^^^^ 

t + 1 A. («„/,)3 aUt+1)' 

Recalling that = 1/ (p„?ii) ~ this implies estimate (15.21) for all t large enough. □ 
Theorem 5.3. Uniformly in t G [0, oo] we have 

lim = u,{t), Um nr'^^feW] = 9,{u,{t)). (5.4) 

n— >oo n— >oo 

/n particular, for t = oo this yields 

Um nr'i?2(6) = 1, lim nr'i?2(6) = c^- (5-5) 

n— >oo n— >oo 

Proof. Similarly to representation (14.171) . one can show that 



oo 



%'^J[6W]=Pn J] m/i(m)^ J2 alx^e-'^^-^'^'^<'^^\ (5.6) 



k,m=l a;i=l 0<X2<tei 



Assuming that to < t < ti and passing to the limit similarly as in the proof of Theorem |4.2[ 
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we obtain, using (14.121) and making the substitution X2 = sxi, that 



k,m=l 

oo ^ oo 



2C(3) p(*)d^7;(n) 

' — (5.7) 



C(2)«3 7o g''{u) 



Similarly, 



oo „ „ 

lim n^^E2[^2it)]= V m/i(m) // e-'=™<^''^W^)Mxidx2 



Erri^im)— — -- / 
k,m=l ^ I -Jta 



2C(3) 



«7(*) q 

ds 



C(2) ^ «:3^;'(n,(s)) 



<(«) 



Finally, the uniform convergence in (15.41) can be proved similarly as in Theorem 15. II □ 
For the future applications, we need to estimate the rate of convergence in (15.51) with suffi- 
cient accuracy. To this end, we require some more smoothness of the function g^^. 

Assumption 5.1. In addition to Assumptions 12. 1 1 and 14. 1 [ we now suppose that G C7^[0, 1]. 
Theorem 5.4. Under Assumption [STl E]{^j) — rij = 0{nl^^) as n ^ oo (j = 1, 2). 
Proof. Consider (the case ^2 is handled similarly). From (15.61) with t = 00 we have 



k,m=l 



00 00 



where 

^i(^) fiihx,,hx2), /i(xi,X2) := xie-^^'''^^^^))^ (5.8) 

xi=l X2=0 

Repeating the calculations as in (15.71) . we note that 



/ / fi{hxi, hx2) dxida;2 
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so that 

/i(m) 



, ™ 1 ^TL \ J i/ 



k,m=l 



E 



h=ankm K,m—L 



= rii. 

n 

(5.9) 



Hence, we obtain the representation 



E:U-ni = ^^Mnkm), (5.10) 



fe,m=l 



where 



Ai(/i) := F^{h) -II /ixs) dxidx2. 

J Jm.\ 

Using that 5i(t) > 5* > and p„ < 1/2, we have 



(I — Q-h5^/2)3 
xi = lx2=0 ^ ' 



Hence, Fx{K) = 0{h~^) as h ^ and Fi{h) = 0{h^^) for any /3 > as /i ^ +00. 
Therefore, the function Fi{h) is well defined for all /i > and its Mellin transform 



00 



Mi(s) := / h'~'Fi{h)dh (5.11) 

(see, e.g., [[HI Ch. VI, §9]) is a regular function for > 2. From a two-dimensional version 
of the Miintz formula (see [|7l Lemma 5.1]), it follows that Mi(s) is meromorphic in the half- 
plane > 1 and has the single (simple) pole at point s = 2. Moreover, for all 1 < < 2 

POO 

Mi(s) = / h'^^Ai{h)dh. (5.12) 
Jo 

The inversion formula for the Mellin transform [[T6l Theorem 9a, pp. 246-247] yields 



1 



c+ioo 



Ai(/i) = — / h~'Mi{s)ds, l<c<2. (5.13) 

In order to make use of formula (15.131) . we need to find explicitly the analytic continuation 
of function (15.1 II) to the strip 1 < 3fts < 2. Let us use the Euler-Maclaurin summation formula 

°° POO 1 POO 

5^/(x)=/ /(x)dx + -/(0)+ / B,ix)f'{x)dx, (5.14) 

x=0 ^ -^0 



where -Bi(x) := x — [x] — 1/2 and [x] is the integer part of x. In view of Assumption 15. II 
and equations (12.111) . (14.211) . we can apply this formula to the sum over X2 in (|5.8I) . Using the 
substitution X2 = txi/pn, we obtain 



foo 1 °° -const-/i 

F,{h) = J2hxi e-'^<^''^(^(^))> dx2 + ^ 5^ /ixie'^^'^'^i^o) + 0(1) - ^ 

xi = l Xl=l 



^ 00 consth. 



xi=l 



I e-^^^'^«dt + 0(l)^-^, (5.15) 
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where (cf. (02l) ) 

m := + tS^it) = Kg'!^{u,{t)f'\ (5.16) 

Keeping track of only the main term in (15.151) and writing dots for functions that are regular 
for 3f?s > 1, the Mellin transform of Fi{h) can be represented as follows 

1 /"OO / °^ />oo \ 



^ oo 



pn ^^^^ Jo \J0 



oo / POO 



h' e-''"i'^(*) dh]dt + - 



Pn X{ Jo Hi) 



\s+l 



dt + 



where 



-C(s-l)r(s + l)^(s) + --- , (5.17) 

Pn 



POO ^ 



Recalling (12.91) . function (15.161) may be rewritten in the form 



llj[t) = KK^itf-' Vl + t^, tQ<t< ti, 



and Assumption 12. 1 1 implies that the function \l/(s) is regular if 3f?s > 0. Furthermore, it is 
well known that the gamma function T(s) is analytic for > ffT3l §4.41, p. 148], whereas 
the zeta function ({s) has a single pole at point s = 1 [[El §4.43, p. 152]. It follows that the 
right-hand side of (|5.17l) is regular in the strip 1 < JJs < 2 and hence provides the required 
analytic continuation of the function Mi{s) originally defined by (15.1 II) . 

Setting h = ankm and returning to formulas (15.101) and (15.131) . we get for 1 < c < 2 

\ r-^MAsKis+ll^^^ ^^^^^^ 



Using that C(s) ^ for 3?s > 1, we can transform the contour of integration 3?s = c in (|5.18l) 
to the union of a small semi-circle s = 1 + re^* (— 7r/2 < t < n/2) and two vertical lines, 
s = 1 ± it (t > r). Furthermore, studying resolution (15.171) . one can show that Mi(l ± it) = 
0(|t|"^) as t — oo. As a result, the right-hand side of (15.181) is bounded by 0(a~^). Thus, 
the proof of the theorem for is complete. □ 

6. Asymptotics of higher-order moments 

Throughout this section, we suppose that Assumptions 12.11 and 14. 1 1 hold. 
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6.1. Second-order moments 

Let Kz := Cov(^,^) = be the covariance matrix (with respect to the measure QJ) of the 
random vector ^ = XIxga- ^'^(•'^)- Since the random variables ^{x) are mutually independent, 
we see using (|3.5|) that the elements Kz{i,j) = Cov{^i, ^j) E {1, 2}) of axe given by 

oo 

Kz{h3) = y^^XiXjVai[iy{x)] = ^ x^ajj ^ fcz''^. (6.1) 

xex xex k=i 

Theorem 6.1. n oo, 

K, = 3K-'n/^{l + o{l))B, (6.2) 
where the elements of the matrix B := {Bij) are given by 

Bii= / M/.^' Bi2 = B2i= / -TTT-TTTT' ^22= / —777-7775- • (6-3) 



9';{uY/^' Jo Jo g'^W 

Proof. Let us consider Kz{l, 1) (the other elements of Kz are analyzed in a similar manner). 
Substituting (14.41) into (16.11) . by the Mobius inversion formula (cf. (15.61) ) we obtain 

00 

k=l x^X 

00 00 00 

= J2 km^f^im) J2 Yl yle-'"'''"^y''^<y^\ (6.4) 

k,m=l 2/1=1^2=0 



Arguing as in the proof of Theorems 14. 2l and 15. 3[ we obtain 



cx> 00 



*i ds 



n— >oo 

a;i = l 2:2=0 ^+ 



Returning to (164)) and using (|4.11l) . (14.121) . we get 

'•'^ ds 3 /"^ dw 



4 • 



C(2) Ao 



and the first formula in (16.31) follows, since a„ = (p„ni)~^/'^ and p„ — t- 1 as n — )■ 00. □ 
Lemma 6.2. A5 — 00, 

det A, . f iV f /' ^ /' 2» - f /' ?M^V^ (6.5) 



Proof. The proof readily follows from Theorem [671] D 
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From Theorem l6.1l and Lemma [6^ it follows (e.g., using the Cauchy-Schwarz inequality) 
that the matrix is (asymptotically) positive definite; in particular, det > and hence 
is invertible. Let Vz = Kz^^"^ be the (unique) square root of K~^, that is, a symmetric, positive 
definite matrix such that = . Recall that the matrix norm induced by the Euclidean 
vector norm | ■ | is defined by \\A\\ := sup|^|^]^ \xA\. We need some general facts about this 
norm (see [|7l §7.2, pp. 33-34] for simple proofs and bibliographic comments). 

Lemma 6.3. If A is a real matrix then \\A^A\\ = || 

Lemma 6.4. If A = {aif) is a real d x d matrix, then 

^E4<pr<E4- (6-6) 

i,j=l i,i=l 

Lemma 6.5. Let A be a symmetric 2x2 matrix with det A ^ 0. Then 

We can now prove the following estimates for the norms of Kz and Vz. 
Lemma 6.6. As n ^ oo, we have 

\\Kz\\-nf, \\Vz\\^nf/\ (6.8) 
Proof. Using Theorem 16. II and the upper bound in Lemma [641 we get 

\\Kzf < Kz{l, If + 2Kzil, 2) + Kzi2, 2f = 0{v?f). (6.9) 
On the other hand, by Theorem 16. II and the lower bound in Lemma [64] 

WK^W" >]^{Kz{l,lf + Kz{2,2f) 

Combining (16.91 ) and (16.91) we obtain the first estimate in (16.81) . 

Further, Lemma 16.31 implies that 1 1 K | P = 1 1 ^ 1 1 . In turn. Lemma 16.51 yields 1 1 ^ 1 1 = 
II A'^ II / det Kz, and it remains to use Lemmas [6^ and [6^ to obtain the second part of (16.81) . □ 

6.2. Asymptotics of the moment sums 

Denote ^^{x) := v{x) — E2[v{x)\ (x G X), and for g G N set 

m,{x) := E][u{xy], /i,(x) := E]\uoixy\ (6.11) 

(for notational simplicity, we suppress the dependence on 7 and z). 

The following two-sided estimate of fj.q{x) can be easily proved using Newton's binomial 
formula and Lyapunov's inequality (cf. ^ Lemmas 6.2 and 6.6]). 
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Lemma 6.7. For each g G M and all x E X, 

/i2(x)«/2 < < 23^^(a;). (6.12) 

Next, we need a general upper bound for the moments of geometric random variables 
proved in [|71 Lemma 6.3]. 

Lemma 6.8. For each g G N, there exists a constant Cq > such that, for all x E X, 

^<ii^) < TT^^- (6-13) 

Using estimate (16.131) and repeating the calculations in the proof of Lemma 6.4 in [|71, one 
obtains the following asymptotic bound. 

Lemma 6.9. For each g G N, 

^ |x|'?mg(x) = 0(1) n^^+^^/^ n -4 oo. 

Lemma [6^ together with bounds (I6.12|) and Theorem 16. 11 imphes the following asymp- 
totic estimate (cf. [|7l Lemma 6.6]). 

Lemma 6.10. For any integer q > 2, 

\x\'^fiq{x) X n^^'^'^^^^, 77, — oo. 

Using Lemma [6. 101 and the lower bound 5j{t) > 5* (see (|4.6I) ). the next asymptotic bound 
is obtained by a straightforward adaptation of the proof of a similar result in [|7l Lemma 6.7]. 

Lemma 6.11. For each q eN, 

E2\ir-E2{ir)\'' = 0{nl''^^), n ^ oo. 
Finally, let us consider the Lyapunov coefficient 

L,:= IlKf J^klVsla:), (6.14) 

xex 

The next asymptotic estimate is an immediate consequence of Lemmas 16.61 and 16 . 1 0[ 

— 1/3 

Lemma 6.12. As n ^ oo, one has . 



7. Local limit theorem 

The role of a local limit theorem in our approach is to yield the asymptotics of the probability 
Q1{^ = n} = Ql^Un) appearing in the representation of the measure P2 as a conditional 
distribution, P,7(-) = QZ{- |i7„) = Q2(-)/gi(i7„). 

As before, we denote := E2{C) ^^id Kz := Cov(^, ^) = EJ. Let /o,/(-) be the density 
of a standard two-dimensional normal distribution A/'(0, /) (i.e., with zero mean and identity 
covariance matrix). 

Then the density of the normal distribution J\f{az, K^) is given by 

fa^,KM = idetKzr'/^fo,i{ix - az)Vz), x G (7.1) 
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Theorem 7.1 (Local limit theorem). Under Assumptions 12.11 and \A. 1 [ uniformly in m G 

QZ{i = m} = fa.,KAm) + 0{nf^'). (7.2) 

Let us make some preparations for the proof. Recall that the random variables {i^{x)}xex 
are mutually independent and have geometric distribution with parameter z^, respectively. In 
particular, the characteristic function (p,y{t) := _E'J(e'*'^) of u{x) is given by 

1 — 

1 — z^e^^ 

hence, the characteristic function (^^(X) := -^^(e'^'^''^^) of the vector ^ = J2x£X ^^{^) reads 

^,(A) = n ^^((^>A);x) = n Y^^S^y (V-4) 

Let us start with a general absolute estimate for the characteristic function of a centered 
random variable (for a proof, see E Lemma 7.10]). 

Lemma 7.2. Let ^^^{t] x) := E2{e^^'^'-^^^^) be the characteristic function of the random vari- 
able Vo{x) := z/(x) — £'2[i^(x)]. Then 

|(/p,o(t;x)| <exp|-|/i2(x)t2 + i/i3(x)|t|='}, t G R. (7.5) 

The next lemma provides two estimates (proved in [|71 Lemmas 7.11 and 7.12]) for the char- 
acteristic function v^^ol^) •= -E'2'(e'^'^'^°'') of the centered vector '■= ^ — clz = Zlxe-v ^^o(a;). 
Recall that the Lyapunov coefficient is defined in (|6.14l) . and Vz := Kz^^"^ . 

Lemma 7.3. (a) For all \ EM?, 

|<^5o(AK)| < exp{-i|A|2 + \Lz\\\^]. (7.6) 

(b) // |A| < L-^ then 

^^,{XVz)-e-\'\'/'\<lQLz\\\'e-\'\'". {1.1) 

The next global bound is obtained by repeating the proof of Lemma 7. 14 in [|7]|. 
Lemma 7.4. For all A G M^, 

|<^?o(A)| <e-'"(^\ (7.8) 

where 

J„(A) ■■=\Y^ e-^r.(&,x) (^^ _ ^^g^^^ > Q (79) 



We can now proceed to the proof of Theorem 17. 1[ 
Proof of Theorem W7\\ By the Fourier inversion formula, we can write 



Ql{^ = ^} = ^^^e-<^''-'^^>(^5„(A)dA, m G Zl, (7.10) 
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where := {A = (Ai, A2) G : |Ai| < tt, IA2I < vr}. On the other hand, the characteristic 
function corresponding to the normal probability density fa^,K^{x) (see (17.11) ) is given by 

¥^a.,x.(A) = e^<^'-)-l^^^"^l^/^ A6M^ 
so by the Fourier inversion formula 

Note that if [AK."^! < then, according to Lemmas [6.61 and [6 . 1 2[ 

|A| < \XV~'\ ■ ||V;|| < L~'\m = 0{n-'^'') = 0(1), 

which of course implies that A G T^. Using this observation and subtracting (|7.1 II) from 
(17.101) . we get, uniformly in m G , 

\Q:{^ = m}- /a.,x.(m)| <h + l2 + h, (7.12) 



where 



4vr^ J {\:\\v.-^\>L-^} 



47^2 

By the substitution A = yVz, the integral Ji is reduced to 



= 0(l)(detK,)-V2L, / |y|Vl^l'/M2/ = 0(n7^/'), (7.13) 

on account of Lemmas 16.21 [6?T2l and iTJl' b). Similarly, again putting A = yVz and passing to 
the polar coordinates, we get, due to Lemmas [6.21 and [6.121 

/ r e-^'/2 dr = 0(n-'/') e'^^^'/^ = o(n-'/'). (7.14) 



I det Vz 



27r 

Finally, let us turn to /a. Using Lemma [74l we obtain 



/a = 0(1) / e--^"(^)dA, (7.15) 

where Jn{X) is given by (17.91 ). The condition |AV"^^^| > L^^ implies that |A| > ^/2r]an and 
hence max{|Ai|, IA2I} > rjOn, where 77 > is suitable (small enough) constant. Indeed, 
assuming the contrary, from (14.51) and Lemmas [6.6l and [6.12l it would follow 

1 < Lz\\V-'\ < LzTjajKzW'/^ = Oil]) ^0 as r/ ; 0, 
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which is a contradiction. Hence, estimate (17.151) is reduced to 

/g = 0(1) U + I J e--^-^^) dA. (7.16) 

\J\\i\>rjan J\\2\>r]OinJ 

Note that, by Assumption 12.11 and formulas (14.211) . the functions 5i{t), 52{t) are bounded 
above, sup^ 6j{t) < 6* < oo. Hence, ( 17.91 ) implies 

JnW > ^e"""'^*(^'i+^^)(l - cos(A,x)). (7.17) 

To estimate the first integral in (17.161) . by keeping in summation (17.171) only pairs of the 
form X = (xi, 1), xi e Z+, we obtain 

JniX) e"""^* > V e"""-^*"^ (1-3? e'(^^"^+^^)) = 3ft — 

xi=0 ^ ^ 

^ ^ \ in 1 Q\ 

~ \ — Q-OLn |X — g-an+iAi| ' 

because 'Ru <\ii\ for any u G C. Since ?7«,j < |Ai| < vr, we have 

Substituting this estimate into (17.181) . we conclude that Jri(A) is asymptotically bounded from 
below by C(?7)a~^ x n^'^ (with some constant Cij]) > 0), uniformly in A such that rjan < 
|Ai| < vr. Thus, the first integral in (17.161) is bounded bv 0(1) exp(— const • nY^) = o(n^^^^). 

Similarly, the second integral in (17.161) is estimated by reducing the summation in (17.91 ) to 
that over x = (1, X2) only. As a result. Is = o{n^^^^). Substituting this estimate, together with 
(17.131) and (17.141) . into (17.121) we get (ITU) , and so the theorem is proved. □ 



Corollary 7.5. In addition to the conditions of Theorem 17. 1[ suppose that Assumption 15.11 
holds. Then, as n ^ 00, 

Q]{e = n} X n^'/l (7.19) 

Proof. By Theorem 15. 4[ = E]{^) = n + 0{nl^^). Together with Lemma [6^ this implies 

\{n — az)Vz\ < \n — az\ ■ \\Vz\\ =0(1). Hence, by Lemma [6^ we get 



ZTT 



and (17.191) now readily follows from (17.2b . □ 



8. Limit shape 

Throughout this section we work under Assumptions 12.11 14.11 and 15.11 Let us first establish 
that a given curve 7 G is indeed the limit shape of polygonal lines F E U with respect to 
the measure QJ (under the scaling F n^^F). 
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Theorem 8.1. For any 5 > 0, 

lim QUren : dcin^^r.-i) <e] = l. 

n— >-oo 

Proof. In view of Theorem l5.1[ we only need to check that for each e > 

Urn Ql\— sup \lrit) - E1\lr{t)]\ > e\ = (8.1) 



n->oo m 0<t<oo 

Note that the random process 

£0 (t) := Irit) - El[lr{t)] (0 < t < oo) (8.2) 

has independent increments and zero mean, hence it is a martingale with respect to the filtration 
J^t '■= <^{^{x), X G X{t), t e [0,oo]}. From the definition of ir{t) (see (14.21) ). it is also 
clear that ip{t) is cadlag (i.e, its paths are everywhere right-continuous and have left limits). 
Therefore, Kolmogorov-Doob's submartingale inequality (see, e.g., ifTTl Corollary 2.1]) gives 



Q: { sup \iUt)\ > riie] < sup Var[£r(t)] < ^Var(£r). (8.3) 

.o<t<oo J {nie)'' o<i<oo n^e'' 



Furthermore, using decomposition (14.21) and Theorem 16. 11 we have 

Var(£/-) = |x| Var[z/(x)] < ^^(a^i + X2) Var[z/(x)] 

= Var(ei) + Var(e2) = 0{n/^). (8.4) 

Finally, substituting (18.41) into (18.31) . we see that the probability on the left-hand side is bounded 
by 0(r27^/^) 0, which proves (O)- □ 
Let us now prove a limit shape result under the measure (cf. Theorem ll.il) . 

Theorem 8.2. For any e > 

lim {r G i7„ :, ddn^T, 7) < 4 = 1- (8.5) 
Proof. Similarly to the proof of Theorem l8.1[ it suffices to show that for each e > 

7 / i^-i/?o. 



lim P2 < sup \n^^ir{t)\ > e } = 0, 
where the random process ir{t) is defined in (|8.2I) . Recalling formula (|3.6I) . we obtain 



. q:\ sup |£°(t)i>£ni 

sup \£Ut)\ >sm\< ^-^^^ L, (8.6) 

To estimate the numerator in (18.61) . similarly to the proof of Theorem l8.1l we use Kolmogorov- 
Doob's submartingale inequality, but now with the sixth order central moment. Combining this 
with Lemma [6.111 (with g = 3), we obtain 

Q: I sup \fr{t)\ > riA < -^E]\ir - EUir)f = 0{nf). (8.7) 



On the other hand, by Corollary 17.51 the denominator in (18.61) decays no faster than at order 

^-4/3 Together with (18.71) . this implies that the right-hand side of (|8.6I) admits an asymptotic 

-2/3 



hound 0{n^ ) — 0. Hence. Theorem 1 8. 2] is proved. □ 
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